Effects of Non-local Stress on the Determination of Shear Banding Flow 
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We analyze the steady planar shear flow of the modified Johnson-Segalman model, which has an 
added non-local term. We find that the new term allows for unambiguous selection of the stress at 
which two "phases" coexist, in contrast to the original model. For general differential constitutive 
models we show the singular nature of stress selection in terms of a saddle connection between fixed 
points in the equivalent dynamical system. The result means that stress selection is unique under 
most conditions for space non-local models. Finally, illustrated by simple models, we show that 
stress selection generally depends on the form of the non-local terms (weak universality). 



Introduction — Oldroyd's proposal [|| that any sensi- 
ble rheological constitutive equation for a general fluid 
should obey the "admissible conditions" has had a great 
influence on the rheology community. In his "admissi- 
ble conditions" , apart from requiring covariance, Oldroyd 
further imposed the "principle of frame indifference" and 
"locality" . de Gennes has shown that the former is no 
longer true if inertia effects are important Q|. Below 
we discuss how one must extend locality to model shear 
banding flow, which appears in some surfactant and 
polymer solutions. 

Shear banding has long been mooted by the polymer 
rheology community in connection with the spurt effect 
in the processing of linear polymer melts (see, e.g. Ref. 
1^ for a review). It is, however, in surfactant worm- like 
micelle solutions that the phenomena has been firmly 
established |§|,||jJJ^, particularly convincingly through 
magnetic imaging experiments For entangled poly- 
mers or surfactant micelles, as the shear rate increases 
the polymers/micelles gradually align and the fluid shear 
thins. According to reptation theory ||l^,|ll|, the fluid 
shear thins so heavily that a maximum appears in the 
shear stress-shear rate relation at a shear rate of approx- 
imately the inverse reptation time. At high enough shear 
rate, the shear stress is presumed to increase again, re- 
sumably due to local {e.g. solvent) dissipation mecha- 
nisms The steady shear rate curve is qualitatively 
like Fig. 1. For intermediate shear rates where the slope is 
negative, homogeneous flow is mechanically unstable. It 
is found in real systems that, when the controlled shear 
rate increases up to the order of the inverse reptation 
time, the fluid becomes a composite flow profile (shear 
bands), where two or more bands, with alternating low 
and high shear rates, coexist with a common shear stress. 

According to Fig. 1, there is a range of shear stress 
[o'mim o'max] within which the constitutive equation is 
multi- valued (three shear rates for a given stress), so 
that shear bands seem possible for any stress within 
[cmin, O'max]- Howevcr, real systems select a well defined 
stress We therefore need a mathematical condi- 

tion, a selection criterion, which determines the coexist- 
ing stress. There are scattered results showing generi- 
cally that steady state solutions behave very differently 
depending on whether or not the constitutive equations 



are local in space For local models, the final steady 
stress depends on the flow history |l^. Therefore an 
additional selection criterion must be imposed to select 
stress, be it variational [l^jl^ or not pO| , pl| . On the 
other hand, steady state analysis of several non-local 
models, either analytically ||l^,|2^ or numerically p3| , p^ , 
yields a well-defined selected stress. It is interesting to 
know whether non-local effects generally lead to a stress 
selection criterion. 
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FIG. 1. Local constitutive relation cr{-y) for the John- 
son-Segalman (JS) model (Eqs. with e = 0.05. 



In this paper, we first present numerical results for 
the modified Johnson-Segalman (JS) model, as a con- 
crete example of stress selection in a non-local model 
that has been well-studied in local form by many work- 
ers 13 1 . Although perhaps not molecularly faithful 
to any fluid, the JS model is a covariant model which 
posesses a non-monotonic flow curve, providing the nec- 
essary ingredient for a qualitative study of shear banding. 
We then present the main result which shows the general 
link between non-locality and stress selection. Finally we 
use simple examples to illustrate the important fact that 
stress selection depends on the details of the interface 
structure, i.e. it has weak universality. 



Numerical Results for the ModiGed JS Model with 
Stress Diffusion — This model is given by 
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p(af +v • V)v = V • cr (1) 

o- = -pl + 2r/K + S (2) 

S = 2G/« - S/t + DV^S (3) 

where the strain rate is k = (Vv + Vv-^)/2. The first 
equation is the momentum equation, with p, v, cr the fluid 
density, velocity, and stress tensor, respectively. The 
stress comprises the pressure p (I denotes the identity 
matrix), the elastic stress S, and a Newtonian viscous 
stress with viscosity rj. We assume an incompressible 
fluid, V • V = 0, enforced by the pressure. The non- 
Newtonian stress S is governed by Eq. (||), where the 
stress is induced by strain rate k with strength described 
by a modulus G, and r is the stress relaxation time. The 
Gordon-Schowalter convected time derivative is 

S = (a* + V • V) S - Vv"^ • S - S • Vv (4) 
+ {1 - a) CE ■ K + K ■ T,) , 

and describes the "slip" between the polymer and fluid 
extension, with a g [—1,1]. 

Eqs. differ from the standard JS model ^ due 

to the additional stress diffusion term I?V^S in Eq. 
El-Kareh and Leal ||2^ have analyzed the microscopic 
dumbbell model {i.e. the case a = 1) and derived this 
term, representing the diffusion of the stress carrying el- 
ement. Note the important difference that the standard 
JS model is a local model, whereas Eq. (jsj) is not. 

Given a shear stress value a^y, we look for the steady 
planar shear flow numerically. For v — v{y,t)ex, the 
modified JS model reduces to |Q 

a ^ S + ej (5) 
TdtS ^ IDdlS - S -^N + 'j (6) 
TdtN = Dd^N ~N + ^S (7) 

The rescaled variables are defined as cr = y/l — axy/G, 
S = ^xy /G, N ^ [(1 - a)S,, - (1 + a)^yy]/2G, 

D — tD, 7 = 2\/l — d?- TKxy, and e = t]/Gt. Because 
the steady planar shear flow solution of Eq. (||) has zero 
Reynolds number (and most macromolecular dynamics 
are in this limit) we can use the low Reynolds number 
limit {p ^ 0), Eq. (||), to obtain the same steady solu- 
tion. Eqs. (Q]^) are the non-trivial dynamics of S (other 
components are uncoupled). Note that there is only one 
non-dimensional parameter e (the viscosity ratio) in the 
steady state problem. 

Given a, we integrate Eqs. (^0) over time t to see 
whether a steady state banding solution can be reached. 
During the integration, the values of S and A'' at the two 
ends y = ztyo are fixed to those of the high and low strain 
rate branches respectively. At i = 0, the functions S{y) 
and N{y) are chosen arbitrarily as Ci tanhy + C2, with 
the constants Ci, C2 fixed by the boundary conditions. 
The value yo is chosen to be sufficiently large that its 



exact value is irrelevant. We find that a steady (elemen- 
tary) shear band solution can be found only at a specific 
shear stress value, a — ao. If the stress value is too 
large or too small, one of the phases shrinks completely. 
The values cfq as a function of the model parameter e are 
shown in Fig. 2. Note that in an infinite system (Tq does 
not depend on the diffusivity D, which just sets the inter- 
face width. However, anisotropic diffusion (a fourth rank 
tensor D in Eq. may lead to diffusivity dependent 
selection. 
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FIG. 2. The selected stress uq of the modified John- 
son-Segalman model as a function of the viscosity ratio e. 
The limits a,nin{e) and amaxie) are also shown. 

Nontransverse Saddle Connection — We have shown 
above that a non-local model can select the stress sharply. 
Now we show that, in general, a non-local model in pla- 
nar shear flow selects stress sharply (could be uniquely), 
provided that the model is a differential equation and 
satisfies rotational and Galilean invariance. 

1. To begin, we observe that the steady state equa- 
tions, like Eqs. without time derivatives, 
comprise a set of ordinary differential equations 
(ODE) for the independent dynamical variables 
{tpi} {S and N for the JS model), rather than 
partial differential equations (PDE), since only 
differentiation in the velocity gradient direction, 
say y, is present in planar shear flow. The 
phase space for the equivalent first order ODE is 
[{ipi} , {dyipi} , {dyipi} ,...]. Higher order gradi- 
ents in the original PDE would entail a ODE phase 
space with large dimension. 

2. To each homogeneous steady flow (a phase) there 
corresponds a fixed point ([{^*} , 0, 0, . . . ]) in the 
ODE phase space since, by definition, homogeneous 
fiow means the variables do not change with y [p7| . 

3. If the interface width is small compared to the other 
length scales in the problem, e. g. the distance be- 
tween two neighbouring interfaces, it is sufficient 
to consider an elementary shear banding solution, 
which describes, in a planar geometry, a composite 
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flow with a smooth interface separating a single re- 
gion of material in the high shear rate phase from a 
single region of material in the low shear rate phase. 
Mathematically, the elementary shear band is a so- 
lution of the ODE which asymptotically connects 
the high and low shear rate phase fixed points be- 
tween y = ±00. Since the fixed points are both sad- 
dle points (see below) and distinct, the elementary 
shear band solution is a heteroclinic saddle connec- 
tion ||. 

4. Let the attractor and repellor basins |^8| of one 
fixed point, say A, have dimensions a{A) and r[A) 
respectively. Another fixed point B has similarly 
defined a{B) and r{B). Note that if the phase space 
has dimension d (for the modified JS model d = 4), 
then a{A) + r{A) < d and a{B) + r{B) < d. A 
saddle connection joins two fixed points, so that it 
is an intersection of the repellor of one fixed point, 
say A, and the attractor of the other fixed point B. 
We denote its dimension s{A, B) [|o). In the ODE 
phase space, the intersection is at least one dimen- 
sion, so that s{A, B) > 1. There is a trivial inequal- 
ity r{A)+a{B) - s{A, B) < d. If equality holds the 
saddle connection is called transverse: it is formed 
by a robust intersection of two manifolds, and is 
structurally stable against small perturbations of 
the ODE parameters (in the modified JS model, a 
ande). The case r (A) + a {B)- s{A, B) < d is called 
non-transverse, and the corresponding saddle con- 
nection is structurally unstable. It is important to 
recognise that if the elementary shear band solu- 
tion is a non-transverse saddle connection, a small 
change of the shear stress (a parameter in the ODE) 
will remove the saddle connection, which gives a 
very sharp stress selection. In another words, given 
an existing shear band solution, a stress perturba- 
tion is singular if it is of the non-transverse type 
saddle connection. On the other hand, if there is 
a transverse saddle connection shear band solution, 
one can perturb the stress to obtain a (slightly) dif- 
ferent shear band solution. Therefore, there is no 
stress selection for a transverse saddle connection. 
(Stress perturbation becomes a regular perturba- 
tion.) 

5. We now prove that: if the model has rotational and 
Galilean invariance, a shear band solution in planar 
shear flow must be non-transverse. Let us momen- 
tarily assume that there is a transverse saddle con- 
nection approaching, say, flxed point A at y = —L 
to flxed point B at y = L. Galilean invariance 
allows us to choose an inertial frame in which the 
flow velocity Vx{y = ±L) = zLV. Since the bulk are 
rotationally invariant and the boundary conditions 
are symmetric under rotation by 180 degrees, one 
can rotate the solution around the vorticity (z) axis 
by 180 degrees to obtain a different shear band so- 



lution which obeys the same equation and the same 
boundary condition {e.g. Vx{iL) = ±V). This so- 
lution approaches fixed point B at y = — L to fixed 
point A at y = L. Fixed points A and B have both 
attractor and repellor directions, and are thus sad- 
dle points. Now take the (thermodynamic) limit 
L — > 00 so that the above two solutions become 
saddle connections. The two saddle connections are 
related by the symmetry transformations, so are of 
the same type (both transverse). A contradiction 
now appears, since "there are not enough dimen- 
sions to put two transverse saddle connections in 
the phase space". Formally, the transverse con- 
dition leads to r{A) -\- a{B) ^ d + s{A, B) > d 
and r{B) -f a{A) = d -\- s[b, A) > d, therefore 
a{A) -\- r{A) -f a{B) -f r{B) > 2d. However, from 
a{A) + r{A) < d and a{B) -\- r{B) < d we have 
a(A) -\- r{A) +a{B)+ r{B) < 2d. So the shear band 
cannot be a transverse saddle connection. 

Weak Universality — One important aspect of stress se- 
lection can be illustrated by a simple model inspired by 
the macromolecular character of many systems exhibit- 
ing non-linear rheology. Let the shear stress a be 



(8) 



where the first, non-linear, part associated with macro- 
molecules is sensitive to a locally averaged shear rate 7 
as opposed to the second, Newtonian, part attributed to 
solvent, which is sensitive to the true local shear rate 7;oc- 
We approximate the local averaging from 'jioc to 7 by a 



gradient expansion, 7 
hence 

7ioc — 



[l-R' (7)V2)7 



7io 



and 



(9) 



We anticipate sensitivity of the non-local scale -R(7) to 7 
through distortion of the macromolecular shape. 

In terms of the locally averaged shear rate 7 we have 



(10) 



where f{j) = ^(7) -I- 777 is, as before, the steady (and 
constant) shear rate flow curve as exemplified in Fig.l. 
Multiplying Eq. (^o|) by i?~^(7)Vy7 and integrating y 
across the shear band leads to 
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7(00) ^(^) 

7(-oo) 



-dj ■ 



(11) 



where 7(1/ — ±00) depends on a through f{'j{y = 
±cxd)) = a. Since ^yjiy — ±00) = 0, an interfacial 
profile must satisfy I'(cr) = 0, which is the condition to 
select the stress. 

According to Eq. (|l^), different functions R{j) give dif- 
ferent I'((t), and hence different selection criteria V — 01 
Therefore, two models with different R{j) but the same 
behaviour in homogeneous flow, 7(7), can behave differ- 
ently upon forming shear bands. The simple case i?(7) 
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independent of 7 corresponds to the equal areas construc- 
tion speculated upon by Ref. , and cannot be regarded 
as generic. Stress selection has weak universality, imply- 
ing that impurities or other effects which changes the 
interfacial properties could in principle affect quantities 
like R{j) and hence alter the selected shear stress. 

For equilibrium transitions, with coexistence equations 
analogous to Eq. (|l0|), terms involving gradients would be 
exactly integrable without an integrating factor {1/R^{j) 
in Eq. ^l]) because they arise from a functional derivative 
of a coarse-grained free energy. The resulting interface 
solvability condition (i.e. a Maxwell construction) is in- 
dependent of the detailed gradient terms, in contrast to 
the weak universality discussed above. Discussion — We 
close with a few comments. First, we have not proven 
existence of shear banding solutions. There are physi- 
cal phenomena where coexistence between bulk states is 
prohibited because it is unfavourable to form an inter- 
face, and robust hysteresis can be expected. Gel swelling 
provides a beautiful example. 

Secondly, the non-transversality condition of the sad- 
dle connection is, strictly speaking, weaker than unique- 
ness. There are two possible situations where uniqueness 
may fail. The first situation happens when, due to an ac- 
cidental situation, the attractive basin of one fixed point 
and the repelling basin of the other fixed point move so 
as to maintain the saddle connection upon changing the 
shear stress cr as a control parameter. The second situ- 
ation happens when more than one isolated stresses are 
selected, i.e. uniqueness fails globally. These have negli- 
gible chance to be realised in a model. Should it happen, 
one may ask for a physical argument for the degeneracy 
to be sure that it is not a mathematical artfact. Gen- 
erally, uniqueness of the stress selection can be expected 
for models with gradient terms 

An interesting question is whether the stress selected 
by non-local effects can be obtained by a variational prin- 
ciple. A conventional variational principle, like the one 
used for the equilibrium phase transitions, relies on the 
volume contribution to a universal functional and gives a 
criterion insensitive to the interfacial structure p9[ |, be- 
cause the interface only contributes a vanishing fraction 
if the total volume is assumed to be large. The model 
Eq. (|o|), illustrates that the way non-local effects se- 
lect stress is different from a variational principle based 
on such a universal functional. Therefore the obvious 
choices of either free energy or entropy production can- 
not generally represent the selection criterion posed by 
the non-local effects. 

After this manuscript was submitted, Yuan published 
a similar modification to the JS model pT[ |. 
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